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Abstract 

I describe a method to calculate a class of three-loop selfenergy 
diagrams for arbitrary internal masses and external momentum. This 
method combines analytical results and numerical integration, and is 
suitable for implementation in a computer program to produce fast 
and accurate results. For the class of Feynman diagrams considered 
in this paper this method leads to a two-fold integral which needs to 
be evaluated numerically. Results are given for a range of masses and 
external momentum. 



1 Introduction 



Perturbation theory based on Feynman diagrams continues to be the main 
tool for calculating observable quantities which can be measured in high en- 
ergy experiments. Unfortunately, the difficulty of the calculations increases 
dramatically with the loop order. 

Some electroweak parameters are being measured in LEP/SLC experi- 
ments at the per mille level. To match such accuracy, full two-loop elec- 
troweak corrections may be needed in the future [P. Higher loop contribu- 
tions are also needed when one deals with relatively large couplings. Here 
higher order corrections can be used to increase the accuracy of the cal- 
culation and to estimate the validity range of perturbation theory. In the 
electroweak sector potentially large radiative corrections may be induced for 
instance by the large Yukawa coupling of the top quark, and also if the mass 
of the minimal Higgs boson turns out to be large. One can also imagine 
extensions of the symmetry breaking sector of the standard model where 
strong interactions may play an essential role |2[. 

In calculating radiative corrections, the difficult step is usually the eval- 
uation of the scalar integrals. 

At one-loop order this is in principle a completely solved problem. How- 
ever, in practice the evaluation of one-loop scalar integrals in the general 
mass case can be tricky because one has to make sure the functions involved 
remain on the physical Riemann sheet, and because of the need to control 
potentially large numerical cancellations |Q. 

At two-loop level, well-known results by van der Bij and Veltman guar- 
antee that all vacuum diagrams with arbitrary masses can be expressed 
analytically in terms of Spence functions If the external momenta do 
not vanish, this is in general not possible any longer. In fact, the two- 
loop selfenergy diagram with three propagators is known to be related to 
the Lauricella function, which is a generalization of the hypergeometric se- 
ries [^]. Such functions are not straightforward to evaluate. Typically they 
are expressed by multiple series whose convergence may be poor for certain 
choices of the parameters, for instance near a threshold. Therefore one has 
to rely at the two-loop level at least partly on numerical techniques. A com- 
bination of analytical and numerical methods exists which can be used in 
principle to calculate any two-loop diagram. There exist also a large num- 
ber of techniques which work for certain topologies or for special choices of 
parameters. On the evaluation of two- loop scalar integrals in the general 
massive case see refs. Q — [pl|] . 



1 



While techniques to handle massless three-loop diagrams exist @, the 
massive case is much more difficult. Apart from a formal expression for 
the four propagator selfenergy diagram no general results exist at the 
three-loop level. Considering the previous remarks on the massive two-loop 
diagrams, there is little hope that analytical techniques can prove useful for 
realistic three-loop calculations of physical relevance in the massive case. 

It is the purpose of this letter to provide efficient methods for evaluating 
certain massive three-loop diagrams. Instead of trying to solve the general 
problem, this paper limits itself to calculating the selfenergy diagram of fig. 
1. With little modifications, this method can be used to calculate some 
other diagrams, of the type shown in fig. 2. 

The scalar integral corresponding to the diagram of fig.l is: 



/(mi, 7722,..., mg; A; ) = 



{q^ — m\)[{q + — "ml] (r^ — m1)[{r — k)"^ — m|] 

(1) 



(p2 — mj)[[p + ky — mg] [{p ~ (lY ~ ^V\\iP + '')^ ~ ""^i] 

To calculate this diagram one starts by performing the integrals over q 
and r analytically, and leaves the p integral to be performed numerically 
after an examination of the analytical structure of the integrand. 

The q and r integrals are essentially two one-loop triangle subdiagrams. 
Of course, in our case two of the external momenta of the triangle subdia- 
grams are functions of the loop momentum p of the three-loop diagram. To 
fix the notations, let us denote: 



r^t 2 2 2 2 2 2\ 



d!^p 



(p2 _ m\)[{p -P2Y - ml] [{p + pif 



ml] 



(2) 



Such triangle graphs are expressible in terms of Spence functions for any 
choice of the internal masses and external momenta. Explicit formulae can 
be found for instance in refs. [^] or |18]. These expressions were encoded in 
the program FF |19]. This program will be used in the following to calculate 
the diagram of fig. 1 in the general mass case. However, the rather lengthy 
expression of the function C simplifies considerably for special choices of the 
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masses. Such is for instance the massles case, mi = m2 = rns = 0. In this 
case a formula can be derived which is better suited for numerical evaluation 
than the general case formula used in FF. Equivalent formulae were obtained 
previously by other authors (see for instance ref. fl^), but the expression 
which will be derived in the following has the advantage of remaining on the 
right sheet for the parameter range needed in this three-loop calculation. 
We therefore sketch the derivation of this formula. 

One introduces Feynman parameters, performs the integral over the loop 
momentum, expands in e = n — 4, integrates over one Feynman parameter, 
and obtains: 



C(0,0,0;p?,pi,p^) = -i— / 

P2, -^0 



dx 



1 



X(l-x)(l-/i2 



log ^? 



(3) 



where 



ax + 6(1 — x) 
x{l — x) 



1 ' 



vl 
vV 



(4) 



It suffices to keep only the 0{^^ terms in the expression of C because 
the three-loop diagram of eq. 1 is obviously ultraviolet convergent. The 
integral of eq. 3 can be carried out analytically and gives: 



C7(0,0,0;p?,pipi) 



Sv\ ) + 'Sp( j + Tlog — + -log — 

v\ Ul 4 4 u\ 



1 2 ""1 1 , 2 ^2 , TT^ 

+7 log 7 log — + 

4 Vl 4 f2 



(5) 



where 



Ml,2 

A 



^(l + fe-a± VA) 

^(l-6 + a±\/A) 
1 - 2(a + 6) + (a - hf 



(6) 
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This function is similar to the finite part of the two-loop integral 
^(mi, m2, TTT-s; 0) of ref. or (mi, mi|m2|m3) of ref. Q. This is not 

surprising, since there is a relation between one-loop massless diagrams 
evaluated at finite external momentum and massive two-loop vacuum dia- 
grams. Other relations of this type were obtained for instance in ref. [p!^]. 

One can now write the three-loop diagram of fig. 1 in the following form: 



I{mi,m2,... ,ms;k ) = 

/ ^^^7^2 2T7^2 2T^("^l'"^2>"ii;P2^ 

J {Pf - mj)[P^ - mil 

X C {ml ml ml Pi, Pi k^) , (7) 



where 



Pi = [p + kf . (8) 

This expression has the structure of a one-loop integral with some vertex 
form factors. As already mentioned, one does not need to keep terms of 
order e and higher in the expression of C because / is ultraviolet finite. One 
can choose the external momentum to have vanishing space components, 
k = (v^, 0). Let us denote the time and space components of the loop 
momentum by p = {po^p)- It is obvious that the integrand in eq. 7 is 
independent of the direction of p. Therefore one can readily perform the 
angular integration over the direction of p. The result reads: 



/ 2 2 
(mi, 7722,, 



mg; 

oo 



dpQ / dp 



{Pi 

xC{m\ 



ms 



2 2 

4, m5, 



m2)[P| 

"^6' P2 7P1 1^ 



■C{mlmlmlPi,Pi,k^ 



(9) 



where 



pI 



pl + 2poVl^ + k^ 



(10) 



and p =\P\- 

In this way one has obtained a two-fold integral representation of the 
three-loop diagram shown in fig. 1. 

Some remarks on the numerical integration are now in place. 
The possible singularities of the integrand in the integration domain 
(pO) p) £ (—00, +00) X (0, +00) are given by the solutions of the equations: 

Po — — ruj + irj = 
pl + 2pQk + k'^ - p^ -ml + ir] = , (11) 

and by the thresholds of the vertex functions C . The threshold behaviour of 
the functions C does not pose special problems with respect to the numerical 
integration of expression 9 because the integrand remains finite at these 
points. One also has to keep in mind that for special values of the mass and 
momentum parameters, additional problems may come from the infrared 
singularities of the function C . C {m\,m}2,rn\]p\,P2,p\) displays an infrared 
singularity for p\ = m^, p\ = m\ and m\ = and cyclic permutations. 

The simplest and most efficient way to avoid the singularities associated 
with eqns. 11 is to perform a Wick rotation of the loop momentum p. 
Instead of the usual rotation of the time component pq, it is more convenient 
to rotate the radial space component p with a negative phase e~*"^, with 
a G (0, 7r/2). This way the integrand becomes free of singularities and can 
be integrated easily. At the same time this procedure allows a useful check 
on the calculation because the result must be independent of the actual 
value of a. 

One can also calculate this integral without introducing a complex p. 
To do this, the integral has to be split along the solutions of eqns. 11. 
Without this factorization of the singularities, the adaptative integration 
algorithms are rather inefficient because the singularities are smeared across 
both integration variables po and p. 

As a check on the calculation, one notices that modifying formula 9 into: 

7/ 2 2 2 I 2\ 

/oo roo p2 

-00 Jo {P',-ml)[P--ml] ^'^'^^ 

gives the two-loop selfenergy diagram shown in fig. 3, which was calculated 
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by a number of authors. For nii = m2 = m.3 = 0, 7714 = ms = \k\ = 1 the 
result is: 

J(0,0,0,1,1;1) = 

vr^ [(1.8472631 ± 1.8 • 10"^) +i (3.4451413 ± 1.8 • lO"*^)] , (13) 

in agreement with aheady existing results . The case mi = m2 = = 
rui = 1715 = = 1 gives: 

J(l,l,l,l,l;l) =7r^[ (.923663 ±2.6- 10"^)] , (14) 

which again agrees with the known value of this diagram ^ . 

Turning now to the three-loop case, fig. 4 shows the results for the 
selfenergy diagram of fig. 1 for mi = m2 = = 7714 = ms = mg = 0, 
mj = nis = 1 as a function of the external momentum k'^. It was checked 
that both eq. 5 and the general formula encoded in FF lead to the same 
results. Fig. 5 shows the same diagram in the all-massive case nii = m2 = 

7713 = 7774 = 7775 = 777-6 = 1/2, 7777 = ™8 = 1- 

Obviously, it is straightforward to modify the formula 9 for calculating 
other three-loop diagrams of the type shown in fig. 2. In a similar way, it is 
also possible to introduce more propagators in eq. 12 for calculating a few 
other two-loop scalar integrals. For example, fig. 6 shows the results for 
diagram 2 a) in the case 777i = 7772 = = = 7775 = rng = 0, 7777 = 1. 

To conclude, a simple method was described which allows one to cal- 
culate efficiently a class of three-loop selfenergy diagrams. It works for 
any masses and momentum combinations in the physical region. As a by- 
product, this yields a new method for calculating certain two-loop Feynman 
diagrams. 
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Figure captions 



Fig.l Three-loop selfenergy diagram in tiie general massive case. 

Fig. 2 Some three-loop diagrams which can be calculated with the same 
methods as the diagram in fig. 1. The dots on internal lines are external 
vertices connected to zero external momenta. 

Fig. 3 A two loop sclfcncrgy diagram which can be calculated with the 
same methods as the three-loop diagram of fig. 1. 

Fig. 4 The real (solid line) and the imaginary (dashed line) parts of the 
scalar three-loop function tt^I, with I defined in eq. 1, as a function of 
the external momentum squared k'^. The masses of the internal lines are 
mi = m2 = ma = m4 = = thq = 0, mr = mg = 1. Note that due to the 
definition of / the absorbtive part of the corresponding diagram is the real 
part of /. 

Fig. 5 Same as fig. 4, but for the all-massive case mi = m2 = ma = 
m4 = m5 = me = 1/2, my = ms = 1. 

Fig. 6 Same as fig. 4, but for the three-loop diagram of fig. 2 a) in the 
case case mi = m2 = ma = m4 = = uiq = 0, my = 1. 
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